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^ ABSTRACT 

We investigate the possibility to use higher order moments of gravitational 

lensed images in the weak lensing analysis. For this purpose we employ spin-2 

component of HOLICs(Higher Order Lensing Characteristics) developed by us. 

We test the weak lensing analysis with spin-2 HOLICs using actual, ground based 

Subaru observations of the massive galaxy cluster Al689(z=0.183). It turns out 

that spin-2 HOLICs of order up to 8 are sufficiently applicable for weak lensing 

0O ' analysis after correcting PSF anisotropy as well as isotropic PSF smearing. 
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.Ifh ! 1. Introduction 

Weak gravitational lensing is now regarded as a powerful tool to reveal the distribution 
of dark matter of clusters of galaxies but also large scale structure. It also useful to measure 
the equation of state parameter of dark energy by observing cosmic shear. 

In the usual treatment of the weak lensing analysis, the quadrupole moment of back- 
ground galaxy images is used to quantify the image ellipticity. In recent years, there have 
been theoretical efforts to include the next order distortion effects as well as the usual 
quadrupole distortion effect in the weak lensing analysis (Goldberg & Natarajan 2002; Gold- 
berg & Bacon 2005; Bacon et al. 2006; Irwin & Shmakova 2006; Goldberg & Leonard 2007). 
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We have proposed to use some combinations of octopole/higher multipole moments of back- 
ground images with definite spin properties which we call the Higher Order Lensing Image's 
Characteristics (HOLICs), and have shown that there is simple relations between HOLICs 
with spin-1 (spin-3) and the third-order weak lensing effects, or gravitational flexion with 
spin-I (spin-3) (Goldberg & Bacon 2005), and thus HOLICs serve as a direct measure for the 
flexion(Okura, Umetsu & Futamase 2007, In the following we refer this paper as OUF2007). 
We have also developed realistic method for measuring flexion by the HOLICs approach by 
fully extended the KSB formalism to take into account higher-order PSF anisotropy as well 
as isotropic PSF smearing (Okura, Umetsu & Futamase 2008. In the following we refer this 
paper as OUF2008). The method has been applied to actual, ground-based Subaru observa- 
tion of A 1689 and has able to obtain a bimodal feature in the central region of the cluster 
which is not seen in the usual quadrupole weak lensing analysis. 

In this paper we consider a possibility to apply HOLICs method to measure shear, 
particularly cosmic shear. Weak lensing shear analysis has suffered from various sort of noise. 
Here we concern the intrinsic noise which is nothing to do with observational condition and 
method of shape measurement. The intrinsic noise can be reduced by averaging a sufficient 
number of background images. Since cosmic shear induces at most of the order of 1% 
distortion, it is absolutely necessary to reduce any source of noises. Thus not only large 
number density of background images but also wide area of survey is necessary for the 
cosmic shear analysis which is not always possible in the actual observation. Therefore it 
would be very useful to have a method to reduce the intrinsic noise. We will try to do 
this by using information of high order, non-dimension and spin-2 shapes which are called 
as "spin-2 HOLICs". The shear distortion has no dimension and is spin-2, thus the shear 
affects strongly not only complex ellipticity but also spin-2 HOLICs. From the differences 
of inner and outer region of image, spin-2 HOLICs of different order have different intrinsic 
noises at least partially. In this way the intrinsic noise can be reduced by averaging spin-2 
HOLICs of different orders. In this paper, we define spin-2 HOLICs of various order, and 
calculate the reduced shear and PSF correction. Then we evaluate the relative precision of 
spin-2 HOLICs of different order by testing them using STEP 1 and find that Spin-2 HOLICs 
up to order 8 (which is a combination of 2 8 multipole moments of shape) are useful in the 
shear analysis. Finally we apply the method to Abell 1689 cluster. 



2. Basis of Weak Lensing 

The gravitational deflection of light rays can be described by the lens equation, 

p = o-ViKO), 



(i) 
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where 6 and (3 are the angular positions of the image and source, respectively, and ip(0) 
is the effective lensing potential, which is defined by the two-dimensional Poisson equation 
as V 2 ip(0) = 2k(0), with the lensing convergence. Here the convergence k = £ m E~; t is 
the dimensionless surface mass density projected on the sky, normalized with respect to the 
critical surface mass density of gravitational lensing S cr i t = (c 2 D s ) / (AitGDdDds), where Dd, 
D s , and Dd s are the angular diameter distances from the observer to the deflector, from the 
observer to the source, and from the deflector to the source, respectively. By introducing 
the complex gradient operator, d = d± + id 2 that transforms as a vector, d' = de^, with 
being the angle of rotation, the lensing convergence k is expressed as 

k = ±dd*^, (2) 

where * denotes the complex conjugate. Similarly, the complex gravitational shear of spin-2 
is defined as 

7 = 7i + 272 = -ddip. (3) 
Note that a quantity is said to have spin-s if it has the same value after rotation by 2n/s. 



3. HOLICs 

HOLICs are introduced by us(OUF2007) as particular combinations with a definite spin 
of multipole moments of images. Here we review briefly the HOLICs formalism. It is very use- 
ful to use complex moment instead of usual moment of images to define HOLICs(OUF2008). 

First, we define complex displacement as 

X{ = d6 1 + id6 2 (4) 
XZ = (d9 1 +id9 2 Y N+M V 2 (d9 1 - i d9 2 Y N ~ M V 2 , (5) 

where d9i means displacement vector from image center and iV and M means order of d9 
and spin number, respectively. Similarly, we define as complex displacement of source 
image defined by the displacement vector (d/3i) from the source center. 

Complex moments of images having brightness distribution 1(0) are defined as follows. 

z» = J d 2 exumi (6) 

where q[I(0)} is an appropriate weight function and is taken as the brightness distribution 
1(6) itself in this section for simplicity. In this notation, the spin-2 HOLICs of order N is 
defined as /Zq (N is an even number). For example, N=2 is complex ellipticity x? and 
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N=4 is i] and N=6 is t>jj according to the notation in OUF2008. Similarly, we define Z 1 ^}^ 
as complex moments of the source image using d(3i and Yj%. 

In the shear dominant field, the local expanded the lens equation is expressed using the 
complex moments as 

Y' = (l-K)(X±- g Xl*), (7) 

where asterisk * means complex conjugate and g is reduced shear 7/(1 — k). Using the above 
form of the lens equation we have the following simple relation between the intrinsic and the 
observed spin-2 HOLICs. 



(s)N ry N 

V (8) 



Z^ Z» N + 2 



z (s)N Z N . 

If we assume the average of intrinsic spin-2 HOLICs to vanish, the reduced shear is 
obtained directly by observing spin-2 HOLICs as 

/ 2 Z?\ 

Since the maximum value of absolute spin-2 HOLICs are 1, the maximum value of \g\ is 
2/(2 + N) in the this approximation . 



4. Weighted HOLICs and PSF correction for real observation 

In the actual observation, the data have random noise and images are smeared by 
isotropic PSF and distorted by anisotropic PSF. In order to correct these effects, we introduce 
weighted HOLICs and PSF correction. The general treatment of the PSF correction may be 
found, for example, in Kaiser, Squires, & Broadhurst 1995 and Bartellman & Schneider 2001. 
The technical details of the PSF corrections in measuring HOLICs is found in OUF2008. 

We redefine the weighted HOLICs by using an extra weight function W(Xl/a 2 ) with a 
characteristic scale a as 

Z N M = J <P0I{O)X N M W (%p) (10) 

The weighted HOLICs of source Z$ N are defined similarly in the source plane using d(3i 
and Yj%. Using the lens equation as well as the expansion of the weight function 

w ^ _ w ^pv^ m w ^ _ n^xl w , ^ (u) 
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the reduced shear effects for spin-2 weighted HOLICs becomes 

,(s)N 7 N ( AT i O yN+2'\ ( AT 9 7 N 7 N+2>\ 7 N 

+ £0 I n 7V ~ _ Z 4 , Z 4 ] * = £2_ _ ^iV „_ (7 iV * /-, 2 X 

+ -lr*N ) 9 r> ry N + Ory N \ 9 - ryN ^0 9 ^ 4 9 A LZ ) 



Z ( S )N Z N y 2 ' y y ^ 2 ' a 2 Z,f J y " 

where HOLICs with ' are measured with W'(x) = dW(x)/dx instead of W(x). 
Next, the anisotropic PSF is corrected as 

yNiso yNobs 

£2 £2 p^V _pN * (l o\ 

7Niso y Nobs 1 Xq 1 4 Ag; l-"- ; 



yNiso yl 
^0 ^0 



where 



d jv _ 1 / + 2 ) , (-W + 2) 

P ° = Z»\ 8 Z ° + ^ 2_/ ° + 2^ Z ° ) (14) 

p n _ 1 f (N-2)(N-4) N _ 2 (N-2) N , 1 7A r +2 'A n ,x 

^ = Z*\ 8 Z * + ^ 2_ ^ 4 + 2^ 4 )> (15) 

where the quantities with upper subscript "obs" and "iso" mean the quantities before and 
after making the anisotropic PSF correction, respectively, and HOLICs with " are mea- 
sured with W"(x) = d 2 W(x)/d 2 x instead of W(x). In the actual observation, higher order 
corrections using higher order spin quantities such as P^ v (spin-6) might be necessary for ac- 
curate shape measurement. However one can see in next section that the constructed mass 
distribution is accurate enough even in the approximation using up to P± . 

The Xq is spin-2 anisotropic component of PSF (OUF2008), and obtained from HOLICs 
measured from star images as follows, 

1 yNobs 

Zj 2star QgX 
^1 pN yNobs ■ 
r ^Ostar 

Finally, we obtain the relation between spin-2 HOLICs and reduced shear from each order 
of HOLICs after isotropic PSF correction as follows, 

yNiso / / C N \ \ / f H N \ 

2 _ / riN I ^Ostar \ D N \ „ / riN / ^Ostar \ D N 



siN I ^Ostar \ r>N \ n I I ^Vstar \ pN \ * /-i 7 \ 

yNiso r ~~ \ "I pN J M) 9 ~ I W ~ I pN A I 9 U'J 

^0 V V r 0star / / V V r 0star J J 

We have tested our shear estimation by HOLICs using the ready-made simulation from 
STEP, and have obtained the input value of the shear when we use higher order PSF correc- 
tions (Cq and C4 ) in all order of HOLICs measurement and do not use faint images in the 
higher order HOLICs shear measurement. This can be seen in Figure [1] where the results of 
analysis of STEP1 simulations are shown in each order of HOLICs shear measurement. In 
this simulation we used about 10 5 brighter images with 20.3 <MAG< 23.3 and 8.6 arcmin 
~ 2 . In order to use fainter images for higher order HOLICs shear measurement we will need 
higher order polarization matrices to correct PSF appropriate for higher order HOLICs shear 
measurement which will be presented in the forthcoming paper. 
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5. A1689 analysis 

In this section, we show the demonstration of applying the weak lensing analysis using 
spin-2 HOLICs with order iV = 2,4,6,8 to Subaru imaging observations of the massive 
galaxy cluster Abell 1689 at z — 0.183. Abell 1689 is one of the best studied clusters(e.g., 
Tyson & Fisher 1995; King, Clowe, & Schneider et al; Bardeau et a!.; Broadhurst et al. 
2005a; Broadhurst et al 2005b; Halkola et al. 2006; Leonard et al. 2007; Limousin et al. 
2007; Umetsu, Broadhurst, Takada 2007; Umetsu & Broadhurst 2007), and therefore serves 
as an ideal target for testing new method. 

We use Subaru Suprime-Cam i'-band data which has 30' x 25' field and 0" .202pixel _1 . 
We analyzed 3000 x 3000 pixel(or about 10' x 10') which is center of the field, and the seeing 
is 0".88 by FWHM. 

We used IMCAT(Kaiser et al. 1995) and some scripts (K. Umetsu, private commu- 
nication) for image detection and measure position, gaussian radius "r ff ", half-light 
radius "?V and magunitude of detected objects in some stage of analysis pipeline. 

We used 81 star images which have 2.2 < < 2.5 and 21 <MAG< 22.5 for PSF correc- 
tion, and we used 3366 galaxies(about 33 arcmin -2 ) with 2.5 < < 10 and 22 <MAG< 25.5. 
Figure [2] is the plots of the average of spin-2 HOLICs before and after anisotropic PSF cor- 
rection. 

Figure [3] are the reconstructed mass distributions measured with HOLICs of order N = 
2,4,6,8. The smoothing is 0".15 by gaussian radius and the interval of successive contours 
is An = 0.2 and the lowest contour is k = 0.2. Reconstructed B-modes of each order have 
dispersion la ~ 0.065,0.060,0.062,0.073, respectively. Figure 0] is the average of these four 
reconstructed mass distribution and the interval of successive contours is An = 0.2 and the 
lowest contour is k = 0.2. Figure [5] are the B-mode of Abell 1689 measured with HOLICs of 
order N = 2, 4, 6, 8, the mapping status is same as figure [3J The contours separate k b = 0.2 
and purple is k b = 0. These figures show that the absolute values of B-mode are less than 
0.2 in almost fields. Figure [6] is the average of four B-mode, the contours separate Kb = 0.2 
and purple is Kb = 0. Averaged B-modes has dispersion la ~ 0.059. 

These results clearly show that we can reconstruct A1689 mass distribution using the 
higher order spin-2 HOLICs as well as complex ellipticity \- The dispersion of B-mode is 
reduced by combining B-mode of each order. Compared with the case of using only complex 
ellipticity \i the dispersion is reduced by about 10%. We can even detect a stronger sub 
clump (ft > 0.4) by using higher order spin-2 HOLICs than by using only complex 
ellipticity. We can see the suppression of the noise in the averaged mass distribution. 
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The detail analysis of Abell 1689 (e.g. radial profiles and so on) will be presented in the 
forthcoming paper. 

6. Conclusion 

We examined the possibility to use higher order multipole moments to measure the 
weak shear. For this purpose we use spin-2 HOLICs of order N larger than 2. In order to 
apply the actual observational data, we introduced the weighted HOLICs and investigated 
PSF correction. Then we apply the method to the massive galaxy cluster Al689(z=0.183), 
and found surprisingly that spin-2 HOLICs of order N up to 8 may be able to measure the 
mass distribution with sufficient accuracy. The reconstructed mass distributions from higher 
order spin-2 HOLICs show a 2nd peak in the central region, and the dispersions of B-mode 
measured by these HOLICs are of the order of that measured by the usual complex ellipticity. 
The combined B-mode is reduced by about 10% and thus the accuracy of the combined mass 
distribution is improved by about 10% over that measured only by the complex ellipticity. 

The improvement is obtained by the partial independence of the shear measured by spin- 
2 HOLICs of different orders. The higher the order, the outer the region of the image the 
shear is measured. Thus spin-2 HOLICs of different order measure the shear in the different 
region of the image. In general the direction determined by the measured shear near the 
central region of the image is not necessarily the same with the direction determined by the 
measured shear in the outer part. This is the basic reason why we can improve the accuracy 
of the measurement of mass distribution. Of course the weight function of the different 
order of multipoles are overlapped each other, and thus the shear measured by them are 
only partly independent. We show in Figure [7] the difference between the directions of shear 
determined by spin-2 HOLICs of different orders as a function of the difference between the 
magnitudes of the shear determined by spin-2 HOLICs of different orders. We show this 
figure because the difference is important to reduce intrinsic noise, and we will address this 
point more detail in the forthcoming paper. This Figures are made by analyzing the data 
obtained in the Subaru imaging observation of the blank region ELAIS Nl. We use 35521 
galaxies in 0.28deg 2 m 2x(12.6) 2 7r arcmin 2 , corresponding to a mean surface number density 
of 35.6 arcmin -2 . They clearly suggest that the shear measured by spin-2 HOLICs of larger 
order is more independent of that measured by the usual complex ellipticity. 

In this letter we have shown the usefulness of the higher order spin-2 HOLICs in the 
shear measurement. E This will open a new possibility to improve the accuracy of the weak 
lensing analysis without doing further observations. Although we have tested our method 
to weak lensing analysis for cluster of galaxies, we can also apply the method to the cosmic 



- 8- 



shear. It is very interesting to see how we can improve the measurement of the cosmological 
parameters by the cosmic shear using higher order spin-2 HOLICs. 

For this purpose, we need more detailed analysis of this method, for example the test 
of other PSF data of STEP I and STEP II. These will be given in the forthcoming paper. 

We thank K. Umetsu for useful discussion and providing his scripts for the image de- 
tection , T. Yamada for the observation of Nl and T. Okamura for useful discussion. The 
work is partially supported by Research Fellowships of the Japan Society for the Promotion 
of Science for Young Scientists for YO. This work is also supported in part by a Grants- 
in-Aid for Scientific Research from JSPS (Nos. 18072001, 20540245 for TF) as well as by 
Core-to-Core Program "International Research Network for Dark Energy". 
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Fig. 1. — The results of analysis of STEP1 simulations in each order of HOLICs shear 
measurement. The horizontal axis are the true(input) shear, and the vertical axis are the 
difference of estimated and true shear with la errors. Errors of estimated shear are same 
levels and less than la error in All order. 
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Fig. 2. — The plots are the average of observed and aniso-PSF corrected Spin-2 HOLICs. 
Red, green, blue and purple are N=2,4,6 and 8 respectively, and circles are observed squares 
is aniso-PSF corrected Spin-2 HOLICs. We can see anisotropic PSF of all order of Spin-2 
HOLICs were corrected. 
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Fig. 3. — The reconstructed mass distribution in the central 10' x 10' region of A1689 by 
spin-2 HOLICs of various order up to 8. The contours are spaced in units of 0.2(« 3a), 
therefore purple ,blue and green lines mean k = 0.2,0.4 and 0.6 respectively. North 
is to the top, and East to the left. 




Fig. 4. — The averaged mass distribution in the same region of Figj3] over that measured 
by spin-2 HOLICs of order N up to 8. 
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Fig. 5. — The B-mode in the central 10' x 10' region of A1689 by spin-2 HOLICs of various 
order up to 8. The contours are spaced in units of lcr(~ 0.2), and purple is 0, blue is 0.2 
and so on. 
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Fig. 6. — The averaged mass distribution in the same region of Figj5] over that measured 
by spin-2 HOLICs of order N up to 8. The contours are spaced in units of 1<t(~ 0.2), and 
purple is 0, blue is 0.2 and so on. 
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Fig. 7. — Figures show the differences between spin-2 HOLICs of order iV = 4(left), 8 (right) 
and complex ellipticity x- The horizontal axis is the differences between the absolute value 
of spin-2 HOLICs of order N=4(8) and absolute value of complex ellipticity. The vertical 
axis is the angles determined by spin-2 HOLICs of order N=4(8) relative to the direction 
determined by complex ellipticity. 



